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■ Abstract 

. Methods of dynamical systems analysis are used to show rigorously that the presence of a 

' magnetic field orthogonal to the two commuting Killing vector fields in any spatially homoge- 

I neous Bianchi type VIo vacuum solution to Einstein's equation changes the evolution toward the 

■ singularity from convergent to oscillatory. In particular, it is shown that the a-limit set (for time 
direction that puts the singularity in the past) of any of these magnetic solutions contains at 
least two sequential Kasner points of the BKL sequence and the orbit of the transition solution 
between them. One of the Kasner points in the a-limit set is non-flat, which leads to the result 
that each of these magnetic solutions has a curvature singularity. 

; PACS number: 04.20.Dw 
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^ ! 1 Introduction 

o : 

Q>^ ^ Progress in the characterization of the asymptotic behavior of spatially homogeneous solutions to 

0^ ■ Einstein's equation has been made by the use of methods from dynamical systems analysis. (See [0 

I for an overview.) For example, these methods have been used to show that solutions in various 

spatially homogeneous subclasses are convergent or oscillatory in a given time direction. The vacuum 
j I Bianchi type VIq spacetimes are a class of spatially homogeneous solutions which are known to be 

bi). convergent in both time directions One time direction is singular and the other is not. In these 

^ I spacetimes the three linearly independent Killing vector fields can be chosen so that two of them 

commute with each other. In |^ evidence is presented that provides strong support for the conjecture 
that if there is a magnetic field present which is orthogonal to the two commuting Killing vector fields 
. then the approach to the singularity is no longer convergent but instead oscillatory, and moreover, 

■ ■ ' generically mixmaster. This is corroborated by the numerical study reported in Q. 

Rigorous results concerning mixmaster dynamics have been elusive. Classes of spatially homoge- 
neous solutions believed to be mixmaster are Bianchi type IX, Bianchi type VIII and magnetic Bianchi 
type VIq, with or without a perfect fluid in each case. (See ||^, |j for additional mixmaster classes.) 
Not only has the occurrence of mixmaster dynamics remained a conjecture, but it has not even been 
shown that solutions in the classes thought to be mixmaster are oscillatory. Progress toward such 
a result for vacuum Bianchi type IX and vacuum Bianchi type VIII was made in The magnetic 
Bianchi type VIq spacetimes considered in are considered again in the present paper (without a 
perfect fluid) and, building on the results of by using techniques developed in [Q, it is rigorously 
shown that all of them are oscillatory in the singular time direction. 

The methods used to show this are those of dynamical systems analysis. A discussion of these 
methods can be found in and in references cited therein. In particular, the formulation of the 
magnetic Bianchi VIq spacetimes as a dynamical system is taken from Q , and definitions and results 
from 1^ and Q are used throughout. Three types (convergent, oscillatory, mixmaster) of asymptotic 
behavior are discussed in section 2. A review of results from is given in section 3. The proof of the 
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new result that magnetic Bianchi VIq solutions are oscillatory is presented in section 4. A method 
which might lead to a proof of mixmaster behavior in this class of spacetimes is suggested in section 5. 

Knowledge of the asymptotic behavior of a solution can lead to an answer to the question of 
whether or not an extension through the initial singularity is possible. In section 6 it is shown 
that each of these solutions has a curvature singularity. Therefore no extension is possible. In the 
course of the analysis, several other types of spatially homogeneous solutions to Einstein's equation 
are considered. Which of these have curvature singularities is also determined in section 6. The 
solutions considered here which can be extended past the initial singularity, and therefore do not have 
a curvature singularity, include rotationally symmetric solutions of Bianchi type I with a magnetic 
field orthogonal to the plane of rotational symmetry. An extension of these spacetimes is given in the 
appendix. Section 7 consists of concluding remarks. 

2 Types of asymptotic behavior 

In a spatially homogeneous solution to Einstein's equation there is a foliation by homogeneous spa- 
tial hypersurfaces with time coordinate t constant on each hypersurface. Assume that the foliation 
of homogeneous spatial hypersurfaces is maximally extended. Choose t so that t takes on all real 
values. Define a generalized Kasner exponent, Pi{t), to be an eigenvalue of the extrinsic curvature 
of the homogeneous hypersurface at time t divided by the mean curvature of the homogeneous hy- 
persurface 0]. In a Kasner (vacuum Bianchi type I) spacetime the generalized Kasner exponents are 
exactly the Kasner exponents and they are constant in time. There are other spatially homogeneous 
solutions to Einstein's equation, for example the spatially flat Friedmann-Robertson- Walker models, 
whose generalized Kasner exponents are constant in time, but do not satisfy the Kasner relation, 
p1 + P2 + — 1. A convergent solution (in a particular time direction) is one in which each of the 
three generalized Kasner exponents converges to a limit in that time direction [Q. For example, if a 
solution is convergent in the negative time direction, then Pi{t) exists. A solution may be 

convergent in both time directions, in only one time direction, or in neither. 

When a class of spatially homogeneous solutions is studied as a dynamical system using the 
expansion normalized variables of Wainwright and collaborators then a convergent solution is 
one whose a-limit set (for the negative time direction, or oj-limit set for the positive time direction) 
consists of a single equilibrium point of the system. Such an equilibrium point may be a Kasner point, 
which represents a Kasner solution, or it may represent another solution whose generalized Kasner 
exponents are constant so that, as in the Kasner solutions, the evolution is completely characterized 
by the volume expansion. An effect of the normalization by the volume expansion is that each 
equilibrium point corresponds to an entire spacetime and each one-dimensional orbit corresponds to a 
one-parameter family of conformally related spacetimes ||] . Since the volume expansion is zero at the 
moment of maximum volume in spacetimes of Bianchi type IX, the expansion normalized variables 
represent half of the evolution - either the expanding phase or the collapsing phase, and so the limit 
set of a solution in the direction of maximum expansion gives information about the evolution at the 
moment of maximum expansion. A magnetic Bianchi type VIq spacetime is cither always expanding 
or always collapsing, depending on choice of time direction, so the expansion normalized variables 
cover the entire spatially homogeneous evolution. 

An oscillatory solution is one whose a-limit set (or w-limit set) contains at least two equilibrium 
points of the system. In all known cases, the equilibrium points contained in such a limit set are a 
finite number of isolated Kasner points and the limit set is a heteroclinic cycle. A heteroclinic cycle 
is a set of n distinct equilibrium points, qi, and n orbits of the dynamical system, F^, such that each 
orbit is convergent in both time directions and the w-limit set of F^ consists of the point q,; while if 
i < n the a-limit set of F,; consists of the point q^+i and if i ^ n the a-limit set of F^ consists of the 
point qi. 

Qualitative analysis and the results of numerical simulations give the following picture for the 
evolution toward the singularity in generic solutions belonging to the classes of spacetimes believed 
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to be mixmaster. Let be a nonexceptional solution in one of these classes. A sequence of Kasner 
evolutions, qi, approximates the evolution of S. Given qi, then qi^i is determined by the BKL 
map 1^, |l^, |l^, 12 1 . Moreover, when the evolution of S is not approximately Kasner, during the 
transition from qi to qi+i, the evolution of S is approximated by a spatially homogeneous solution to 
Einstein's equation which is uniquely determined by qi. In the case of Bianchi types VIII and IX this 
is a vacuum Bianchi type II solution. In the case of magnetic Bianchi type VIq some of the transitions 
are approximated by vacuum solutions of Bianchi type II and some are approximated by magnetic 
solutions of Bianchi type I |^ . The solution approximating the evolution during the transition from 
one Kasner to the next converges to qi to the future and to g^+i to the past (for time direction that 
puts the singularity in the past). As the singularity is approached the approximation of the evolution 
of S by such a sequence of Kasners along with the transition solutions determined by the sequence 
improves. However, it is not necessarily expected that the evolution of 5* converges to a single BKL 
sequence of Kasners along with the transition solutions determined by the sequence |Q . It may be 
that the evolution always eventually diverges from one sequence and another one becomes a better 
approximation. 

The description of mixmaster evolution just given can be stated in terms of the a-limit set of 
S. The presence of a perfect fluid does not significantly alter the picture. If there is no perfect 
fluid, then each mixmaster class is a four-dimensional invariant set of a dynamical system and the 
solution, S, is represented by an orbit in the four-dimensional invariant set. The boundary of this 
four-dimensional set is in each case the union of the three-dimensional sets listed below and the 
compact two-dimensional invariant set. A, consisting of the Kasner points along with the orbits of 
all transition solutions between sequential Kasner points of the BKL sequence. The boundary of 
vacuum Bianchi type IX consists of orbits of vacuum Bianchi types VIIq (a three-dimensional set), 
II (a two-dimensional set) and I (the Kasners). In this case A is the union of the latter two sets. 
The boundary of vacuum Bianchi type VIII consists of orbits of vacuum Bianchi types VIIq, VIq (a 
three-dimensional set), II and I. A is again the union of the latter two sets. The boundary of magnetic 
Bianchi type VIq consists of orbits of vacuum Bianchi type VIq, magnetic solutions of Bianchi types II 
(a three-dimensional set) and I (a two-dimensional set) and vacuum Bianchi types II and I. In this 
case A is the union of the latter three sets. It is known that the a-limit set of S must be contained 
in the boundary of the four-dimensional invariant set in which it lies ^, H] . If 5 is magnetic Bianchi 
type VIq then it is known that its a-limit set is non-empty since the orbit is contained in a compact 
set of the dynamical system. If S is Bianchi type VIII or IX then until recently it has not been known 
that its a-limit set is non-empty. (It appears that this is no longer the case, due to a recent result 
by Ringstrom showing that such an a-limit set is indeed non-empty [|l3).) To be consistent with the 
picture of mixmaster dynamics, the a-limit set of S should be non-empty, and it should be contained 
in just part of the boundary, obtained by excluding the three-dimensional sets. That is, it should 
be contained in the set A. This is because in order for the approximation of the evolution of S by 
a BKL sequence of Kasners along with the transition solutions determined by the Kasners in the 
sequence to get better as the singularity is approached, the orbit of S must eventually remain within 
any compact set containing a neighborhood of A. Given any point, p, of the boundary that is not in 
A, such a compact set can be found which does not contain p. Furthermore, if it could be shown that 
the a-limit set of S is contained in the set A, then the results in [|j and in the present paper would 
show that an entire BKL sequence of Kasner points along with the transition solutions between them 
must lie in the a-limit set of S. Whether ^ or a subset of A is the a-limit set of generic solutions in 
the mixmaster classes would still be an open question. 



3 



3 Magnetic Bianchi VIq 



The union of the orbits of magnetic Bianchi VIq solutions is the four-dimensional open invariant set, 
U, defined by fl: 



+ 1]^ +iV^ + -i?^ = 1, (1) 



3 

iV_ > 0, Nl<3Nl, and H > 0. (2) 



This formulation is in terms of the shear tensor, the structure constants and the magnetic field in 
an orthonormal frame which is invariant under the action of the group of isometrics. H is obtained 
from the magnetic field, which is orthogonal to the two commuting killing fields. "E^ and S_ are 
obtained from the shear tensor while -/V+ and iV_ are obtained from the structure constants. These 
quantities have been normalized by the expansion scalar. Equation]^ is the Hamiltonian constraint. 
It is equation (2.19) in |^ with 17 = 0. See Q for further details concerning this construction and its 
representation of all solutions in this class. 

The time direction is chosen so that the singularity is in the past (r decreasing). The evolution 
equations are: 

S'+ = -27V2(i + 5]+) + ^i/2(2-S+), 
3 

S'_ = -{2Nl + -H^)J:^ -2N+N^, 

iv^ = {2e+{i + j:+) + 2j:^_ + ^h^)n+ + 6i:^n^, (3) 



N'_ = (2S+(1 + E+) + 2E2_+ i72^)^__^2E_iV+, 



3 
2' 

H' = -(S+(2- E+) - +iV2)i/. 



These equations can be obtained from equations (2.17) in Q by using their equation (2.16) with 
f2 = 0. Evolution satisfying equations |^ preserves the Hamiltonian constraint. 

There is exactly one equilibrium point in U, the point at which — E_ ~ 0, = 0, 
= I and H = ^. In ||] it is shown that the cj-limit set of each solution in U consists of this 
single equilibrium point. That is, all of these magnetic Bianchi VIo solutions are convergent in the 
non-singular time direction. 

The boundary of U, dU, is made up of various sets which are invariant with respect to equations |[ 
Solutions to equations |^ in these invariant sets represent solutions to Einstein's equation in various 
other spatially homogeneous classes. As mentioned previously, these are vacuum Bianchi VIq {H — 0), 
magnetic Bianchi II {Nl = 3Nl) and I (iV_ = 0), vacuum Bianchi U {H = and Nl = 3Nl) and 
the Kasner points {H = and — 0). Note that the magnetic field is not the most general one 
possible. See ^ for other possibilities. Each Kasner point is an equilibrium point and there are no 
other equilibrium points in dU. Three of the Kasner points represent spacetimes which are isometric 
to a part of Minkowski space, and so are called the flat Kasner spacetimes. These are labeled Ti, 
i G 1,2,3 (see figure ||). 

The orbit structure in the two-dimensional invariant sets of the boundary, vacuum Bianchi II and 
magnetic Bianchi I, has been known for some time and references cited therein]. Each of these 
solutions is convergent in both time directions. These are the transition solutions from one Kasner 
point to the next in the BKL map. There are two copies of vacuum Bianchi II in dU . On one copy 
(which will be denoted vacuum Bianchi 11+) = y/3N- and on the other (which will be denoted 
vacuum Bianchi II~) N+ = — \/3^-- The Kasner circle is the boundary of each of these three two- 
dimensional open sets. The union of these three sets and the Kasner circle is the set A for magnetic 
Bianchi VIq. 
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Table 1: Strictly monotonic functions 





magnetic 
Bianchi 11+ 


magnetic 
Bianchi II~ 


vacuum 
Bianchi ff+ 


vacuum 
Bianchi 11^ 


magnetic 
Bianchi f 


L 


decreasing 


increasing 


decreasing 


increasing 


constant 


W+ 


increasing 


indefinite 


constant 


decreasing 


increasing 


W" 


indefinite 


decreasing 


increasing 


constant 


decreasing 



The orbit structure can be seen by considering the behavior of the following three functions on 
solutions. 

L 



2- 




H 




V3 




1 





W+ = ^\ + (4) 

Calculation of the time derivative of each of these functions, using equations ^ restricted to each set in 
turn, shows that on each set one of these functions is constant on solutions, one is strictly increasing 
on solutions and one is strictly decreasing on solutions (see table |l|) . Thus the orbits are straight lines 
when projected into the E+-I]_ plane, and it can be determined which Kasner point is in the a-limit 
set and which Kasner point is in the w-limit set of each solution. The cj-limit set of a solution in 
magnetic Bianchi I consists of a point of /Ci, and the a-limit set consists of a point of /C2 U Ti U /C3. 
See figure |l| for the definition of /C^ . The cj-limit set of a solution in vacuum Bianchi 11+ consists of 
a point of /C2 and the a-limit set consists of a point of /C3 V^T^2,V^¥i\. The cj-limit set of a solution in 
vacuum Bianchi II~ consists of a point of /C3 and the a-limit set consists of a point of /Ci U T3 U /C2- 

The polarized solutions in J7, in which one of the degrees of freedom of the gravitational field is 
absent and the metric in a coordinate basis can be diagonalized , form a two-dimensional invariant 
subset of J7, defined by E_ — and = 0. In it is shown that the a-limit set of each polarized 
solution in IJ is the heteroclinic cycle made up of the following orbits: the Kasner points T\ and Qi 
(see figure |l] for definition of Q;), the single vacuum polarized Bianchi VIq orbit, Fi, whose a-limit 
set consists of the point Qi and whose w-limit set consists of the point Ti, and the orbit in magnetic 
Bianchi I, whose a-limit set is Ti and whose w-limit set is Q\. Thus the polarized solutions in U 
are oscillatory toward the singularity. They are not mixmaster. Their evolution toward the singularity 
is not approximated by a BKL sequence. r2 is a transition solution between two sequential points in 
the BKL sequence, but Fi is not. If a BKL sequence includes the point Ti then it ends there. Thus, 
a polarized solution in U is an "exceptional solution" in this mixmaster class (see section 2). 

In 1^ it is shown that the a-limit set of a non-polarized solution in IJ must be contained in dU by 
considering the function 

3S]2 + 7V2 

Z= ^ tr (5 

which is strictly decreasing on any non-polarized solution in \J . 



4 Magnetic Bianchi VIq solutions are oscillatory 

That the magnetic Bianchi VIq solutions are oscillatory toward the singularity is shown by the fol- 
lowing theorem. 
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Theorem. The a-limit set of a spatially homogeneous solution to Einstein's equation with Bianchi 
type VIo symmetry and a magnetic field orthogonal to the two commuting Killing vector fields con- 
tains at least two sequential Kasner points of the BKL map and the orbit of the transition solution 
connecting them. 

Proof of Theorem. The proof of the theorem will be given in the rest of this section. The result has 
already been obtained for the polarized solutions in [/ so it remains to obtain the result for the 
non-polarized solutions in U. It has also already been shown that the a-limit set of a non-polarized 
solution in U is contained in dU ^ . Therefore the first part of the proof of the theorem is an analysis 
of the orbit structure in dU. It is shown that each limit set of any solution in dU contains a Kasner 
point. From this it immediately follows that the a-limit set of a non-polarized solution in U contains 
a Kasner point, because a non-empty a-limit set is a union of orbits along with the limit sets of each 
orbit. Therefore if an a-limit set is contained in dU it must contain an entire orbit in dU and it must 
contain the orbit's limit sets. Then, through analysis of the orbit structure in the closure of U, U, in 
the neighborhood of a Kasner point, it is shown that if the a-limit set of a solution in U contains a 
Kasner point, it necessarily contains a second Kasner point. 

Proposition 1. Each solution in the boundary of U is convergent to a Kasner point in both time 
directions. If the solution is not Kasner, then the Kasner point in its a-limit set is not the Kasner 
point in its tj-limit set. 

Proof of Proposition 1. This result is already known for the Kasner points, the vacuum Bianchi II 
solutions and the magnetic Bianchi I solutions that lie in dU . 

Consider the three-dimensional open invariant set, /a, made up of all the magnetic Bianchi 11+ 
[Nj^ = ^/3N-) orbits in dU. The function L is strictly decreasing on solutions in J3 and is 
strictly increasing. The existence of a function that is strictly monotonic on solutions in an invariant 
set implies that the limit sets of the solutions must be contained in the boundary ^, Proposition Al]. 
The boundary of is made up of the union of all the Kasner points, the magnetic Bianchi I orbits 
and the vacuum Bianchi 11+ orbits. L is strictly decreasing on the latter solutions and is strictly 
increasing on magnetic Bianchi I solutions. A function that is strictly monotonic on a solution must 
be constant on each limit set of that solution and a non-empty limit set is the union of orbits. 
Therefore each limit set of a solution in I^, must lie in the Kasner circle. For a given value in the range 
of L or there are at most two Kasner points at which the function takes this value. A limit set 
must be connected. Therefore each limit set of a solution in I3 must consist of a single Kasner point. 
The Kasner point in the a-limit set of a solution is not the the Kasner point in the w-limit set, since 
a function that is strictly monotonic on a solution cannot have the same value on both limit sets of 
the solution. 

A similar argument, using the functions L and W~ (see table Q), shows that a solution in the set 
I2, made up of all the magnetic Bianchi II~ ~ — \/3A^-) orbits in dU, is convergent to a Kasner 
point in either time direction, and that the Kasner point in the a-limit set is not the the Kasner point 
in the w-limit set. 

Now consider the three-dimensional open invariant set, Ii, made up of all the vacuum Bianchi VIq 
orbits in dU. It was shown in |^ that each of these solutions is convergent to a Kasner point in either 
time direction. In the present context, since the boundary of Ii is made up of the Kasner points and 
both sets of vacuum Bianchi II solutions, this result can be obtained with the same reasoning used for 
I3 by considering the function E+ which is strictly decreasing on solutions in /i and is also strictly 
decreasing on both sets of vacuum Bianchi II solutions. And again, the existence of the monotonic 
function implies that the Kasner point in the a-limit set is not the Kasner point in the w-limit set. □ 

It has now been shown that each limit set of any solution in dU contains a Kasner point, from 
which it follows, as discussed above, that the a-limit set of a non-polarized solution in U contains a 
Kasner point. To show that this in turn implies that there is another Kasner point in the a-limit set, 
further analysis of the orbit structure in U will be useful. 
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Figure |l| shows the sign of the eigenvalues of the hnearization of the differential equation about 
each Kasner point At each Kasner point the eigenvector corresponding to the eigenvalue Xjj is 
tangent to magnetic Bianchi 1 while the eigenvectors corresponding to the eigenvalues \± are tangent 
to vacuum Bianchi IV^ . The fourth eigenvalue, Ak, vanishes at each Kasner point and its eigenvector 
is tangent to the Kasner circle. From the number of eigenvalues which vanish, the number which 
are positive and the number which are negative it follows that each flat Kasner point has a three- 
dimensional center manifold and a one-dimensional unstable manifold. Each non-flat Kasner point 
has a one-dimensional center manifold, a two-dimensional unstable manifold and a one-dimensional 
stable manifold. A center manifold is not necessarily unique, while the stable and unstable manifolds 
are unique. If the orbit structure on a center manifold is known, the orbit structure of the nonlinear 
system near a Kasner point can be determined. Let Oc be the intersection of a center manifold of the 
Kasner point with some neighborhood, iV, of the Kasner point in U . Let j , k and I be the dimensions 
of a center manifold, the unstable manifold and the stable manifold, respectively. Let x = (xi, ...,Xj), 
y = ■■■,yk) and z — (zi, ...,2;/). Here xi, ...xj are coordinates on Oc, y € R'' and z e RK For some 
positive e, let 

O = {(xi,...,a;j,7/i,...,yfc,zi,...,z;) : x e Oc,0 < yi < e and < < e}. 

Let x(t) satisfy the nonlinear system restricted to Oc. Let j/(t) — yoe^ and z(r) = zqb^'^. Then 
there exists an e and an N such that the nonlinear system on N is topologically equivalent to this 
system on O Theorem 2.3]. (See also jl^, page 99].) Topological equivalence means that there 
is a homeomorphism, H : N O, which maps orbits in onto orbits in O and preserves the time 
orientation on the orbits. 

Let g be a non-flat Kasner point. Then the Kasner circle is a center manifold of q. Since the Kasner 
circle is made up of equilibrium points, then in this case x(t) as defined above is independent of r. 
The solutions on O are therefore known exactly. Then the topological equivalence of the nonlinear 
system on A^ to this explicitly known system on O implies the following. In A^ the Kasner circle is 
intersected by a foliation of three-dimensional leaves. Any orbit of the nonlinear system restricted to 
A^ lies entirely in one leaf. (However, the intersection of an orbit in U with A^ may have points on more 
than one leaf.) One of the leaves contains g, which is its own a-limit set and its own cj-limit set. In 
this leaf there is a one-dimensional stable manifold of q, tangent to the eigenvector corresponding to 
the negative eigenvalue. The oj-limit set of solutions on this manifold is q. There is a two-dimensional 
unstable manifold of q, tangent to the two eigenvectors corresponding to the positive eigenvalues. 
The a-limit set of solutions on this manifold is q. Now it will be shown that the stable and unstable 
manifolds lie in dU (and therefore not in U). 

The one-dimensional stable manifold of q has already been found. If g G /Ci its stable manifold 
is an orbit of magnetic Bianchi I. If g e JC2 its stable manifold is an orbit of vacuum Bianchi 11^. If 
q € JC3 its stable manifold is an orbit of vacuum Bianchi 11^ . 

Now consider the closure of /i, Ii. In this three-dimensional subsystem the eigenvalues of the 
linearization of the differential equation about each Kasner point are X± and Xk, with the same 
eigenvectors as in the four-dimensional system. If g g /Ci and A is a neighborhood of q in Ii, then 
the nonlinear system restricted to A^ n A^ is topologically equivalent to the system on O with I set to 
zero. Therefore, in A" n A^ the Kasner circle is intersected by a foliation of two-dimensional leaves, 
each of which is a two-dimensional unstable manifold of the Kasner point that it contains. Since Ii 
is an invariant manifold, this implies that the unstable manifold of q lies in Ii. 

Similar reasoning shows that if g G IC2 , its unstable manifold lies in 1 2 , while if g G /C3 its unstable 
manifold lies in I3. 

Now let S* be a non-polarized solution in U and let A be its a-limit set. It was shown above that 
A contains a Kasner point. Consider the case that A contains q. It has just been shown that S is 
not in either the stable or the unstable manifold of q. Moreover, analysis of the orbit structure in O 
shows that no other solutions there are convergent to H{q) which shows that S is not convergent to 
q. But since q is in A, then there is a sequence of points on 5* and in A^, S'(t„), which converges to q, 
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with r„ < Tn-i. Let r„ be the orbit in N on which S'(t„) Hes. Note that r„ need not be in the leaf 
containing q, but the orbit structure in every leaf is the same and, in addition, the sequence of r„'s 
approaches the leaf containing q. Choose d < e and define 

Os {{x,yi,y2,z) : x e Oc,0 < Vi < 5 and Q< z <5}. 

Then each iJ(r„) has two points in the boundary of Og in H{N n U), one preceding H{S{Tn)) in 
time (because of the stable manifold) and one following (because of the unstable manifold, but the 
preimage of this point will have r < t„_i)). This gives two sequences of points, one which converges 
to H{ps), with ps a point on the stable manifold of q, and one which converges to H{pu), with p„ a 
point on the unstable manifold of q. Each point in each sequence is the image of a point of the orbit 
of S. Thus a point of the stable manifold of q and a point of the unstable manifold of q are in A. 
For a more detailed argument, see pages 33 and 34]. The reasoning is taken from (H Proof of 
Theorem 4-3] and the proof of a similar result in |]l6| . Appendix]. From this in turn it follows that the 
entire stable manifold, its a-limit set, an orbit of the unstable manifold and its a;-limit set must all 
lie in A. The a-limit set of the stable manifold of q is the Kasner point which follows q in the BKL 
sequence. The stable manifold of q is the transition solution between them. 

The other possibility is that A contains a flat Kasner point T^. li is a center manifold of Ti in U. 

Proposition 2. The w-limit set of any solution in /; consists of the point Ti. 

Proof of Proposition 2. It was already shown that the w-limit set of any solution in li consists of 
a single Kasner point. The analysis of the orbit structure in the neighborhood of a non-flat Kasner 
point, q, showed that the stable manifold of g is a single orbit in one of the two-dimensional boundary 
sets, and that the w-limit set of no other solution in U (besides q itself) is the single point q. The 
remaining possibilities are the flat Kasner points, Ti. 

The function is strictly increasing on solutions in I3 . Therefore the value of on the w-limit 
set of a solution in I3 must be greater than the infimum of on J3 . So Ti cannot be in the w-limit 
set of a solution in . The function L is strictly decreasing on solutions in I3 . Therefore the value of 
L on the w-limit set of a solution in I3 must be less than the supremum of L on I3. So T2 cannot be 
in the w-limit set of a solution in I3 . But the w-limit set of a solution in J3 consists of a single Kasner 
point. The only remaining possibility is T^, so must be the cj-limit point of all solutions in I^. 

Similar reasoning shows that the tj-limit point of all solutions in I2 is T2 . 

To find out what is the case for Ii, first recall the polarized case. This one-dimensional set is a 
single orbit whose a-limit set is Qi and whose w-limit set is Ti. Now consider the (open and invariant) 
set of non-polarized solutions in Ii . The function Z (equation |^) is bounded from below (Z > 0) and it 
is strictly decreasing on solutions in this set. Therefore it must approach a finite value on the o^-limit 
set of one of these solutions. This is impossible at T2 and at T3. However, at Ti, if (3S^ + Nf) 
fast enough along solutions compared to 3A''^ — iV^ then it is possible for Z to approach a finite value. 
Since the w-limit set must consist of a single Kasner point, and all the other possibilities have been 
ruled out, then this must be the case. □ 

Proposition 3. The a- limit set of any solution in Ii consists of a point in /C^. 

Proof of Proposition 3. Given a point g G /Ci , it has been shown that the a-limit sets of solutions 
on a two-dimensional submanifold of U consist of the point q. It has been ruled out that the a-limit 
set of any other solution in U can consist of the single point q. It has been shown that the a-limit 
set of a solution in Ii consists of a single Kasner point. The point Ti is ruled out by the existence of 
a monotonic function on solutions of 7^, since the w- limit set of solutions of Ii consists of the point 
Ti. The remaining two flat points, Tj, j ^ i can be ruled out in each case by considering Hj{Ii n Nj), 
defined in the next paragraph. □ 
Now consider the case that a flat Kasner point Ti is in A. Let Oc be the intersection of Ii with a 
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neighborhood, Ni C U, of T^. Let xi,X2.,xz be coordinates on Oc- Let 

Oi = {{xi,X2,X3„y) : x e Oc,0 <v<e}. 

Let a:;(T) satisfy the nonhnear system restricted to Op and y(T) = yoe^- Let e and A^i be such that 
the nonhnear system on Ni is topologically equivalent to this system on Oi, with homeomorphism 
Hi : Ni ^ Oi. If X is on the boundary of li then H~^{x, y) is on the boundary of U . If a; G Ii and 
< ?/ < e then H~^(x,y) € U. Since there exists a (small enough) neighborhood of Ti such that 
no solution in li lies entirely within the neighborhood and since all solutions in li converge to Ti, li 
plays the role that the stable manifold played for the case that a non-flat Kasner point is in A. Then 
similar reasoning leads to the result that the unstable manifold of Ti and an orbit other than Ti in 
the closure of the center manifold along with the limit sets of these two solutions must be in A. Qi is 
the Kasner point preceding Ti in the BKL sequence, and the unstable manifold of Ti is the transition 
solution between them. This ends the proof of the Theorem. 

5 From oscillatory to mixmaster 

It was found in the previous section that A contains the following: at least one non-flat Kasner point 
(gfc), the stable manifold of q^, its a-limit set [qk+i), an orbit in the unstable manifold of qk and its 
w-limit set (which is a single Kasner point). If qu+i is a non-flat point, then its stable manifold is also 
in A, and so on. Let qu — > Ik+i stand for the union of q^ and its stable manifold and qk+i- Then 

qk -> qk+i qk+2 qk+3 qk+4 ■■■ (6) 

is in A. This BKL sequence of Kasner points can end at a flat point, can be periodic, or can be 
non-ending and non-periodic ||l^, |l^, |l^ . 

For this class of solutions to fit the standard picture of mixmaster dynamics it must be the case 
that orbits in the unstable manifold of qk that do not belong to A are not in A, for generic solutions 
in U. That is, there should be no points of vacuum Bianchi VIq and no points of magnetic Bianchi II 
in A. A polarized solution in U is exceptional. There is an orbit, Fi (see the end of section 3), in its 
a-limit set that does not belong to A. Fi is the orbit of the vacuum polarized Bianchi VIq solution. 
It will now be shown that no other solutions in U have this type of exceptional behavior. 

Proposition 4- No points in the orbit of any vacuum Bianchi VIq solution are in the a-limit set. A, 
of a non-polarized magnetic Bianchi VIq solution. 

Proof of Proposition 4- The function Z (equation ^ is strictly decreasing on non-polarized solutions 
in U, so Z must be constant on A. Z is also strictly decreasing on non-polarized solutions in vacuum 
Bianchi VIq. Since a-limit sets are the union of orbits then no point of a non-polarized vacuum 
Bianchi VIq orbit can be in A. The polarized vacuum Bianchi VIq orbit is also ruled out since the 
value of Z at points of this orbit equals the infimum of Z on the set of non-polarized solutions in U. 
□ 

This means that if some qk G /Ci is in A then the only orbits in its unstable manifold that can 
be in A are the two that are in vacuum Bianchi II. These are precisely the orbits of the transition 
solutions between qk-i and qk- (Given a non-flat Kasner point q, there are two Kasner points which 
can immediately precede q in the BKL sequence.) 

Two possibilities come to mind for showing that no points of magnetic Bianchi II are in A, for 
generic solutions in U. One possibility is more careful analysis of the orbit structure in U in the 
neighborhood of a Kasner point. Another is the following. If two functions analogous to Z could be 
found, both of which were strictly monotonic on (perhaps generic) solutions in U, one also monotonic 
on (perhaps generic) solutions in magnetic Bianchi 11+ and the other monotonic on (perhaps generic) 
solutions in magnetic Bianchi II~, then the same result just obtained for qk G /Ci would be obtained 
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for qk (z IC2 and for qk G /C3. Only two orbits in the unstable manifold of qk would be left as 
possibilities for being in A. These are the orbits of the transition solutions of the BKL map between 
the two possible qk~i^s and qk- One or the other of these orbits would have to be in A. Therefore a 
sequence of the following form would have to be in the a-limit set of generic magnetic Bianchi VIq 
solutions. 

• • ■ qk-2 qk-i qk qk+i qk+2 ■■■ (7) 

This sequence is either periodic or infinite, since, while the domain of the BKL map is the set of 
all Kasner points except for the three flat Kasner points, the range is the set of all Kasner points. 
Therefore the sequence can end (a point is in the range that is not in the domain), but it cannot 
begin (every point that is in the domain is in the range). 

Showing that no points of magnetic Bianchi II belong to A would, furthermore, imply that if the 
sequence ends at a Ti, that is, if Ti is in A, then so is a non-flat Kasner point in any neighborhood 
of Ti. As of now, this result can only be obtained for Ti. 

Proposition 5. If Ti is in the a-limit set, A, of a non-polarized magnetic Bianchi VIq solution, then 
for any neighborhood, N^, of Ti, no matter how small, there is a non-flat Kasner point q such that 

qe N^nA. 

Proof of Proposition 5. In the proof of the theorem, it was found that if Ti G A then an orbit in Ii 
is also in A. But it has now been shown that no points of vacuum Bianchi VIo are in A. Therefore, 
an orbit in vacuum Bianchi II is in A. Furthermore, the reasoning in the proof of the theorem shows 
that for any neighborhood Ng of Ti that is small enough, there is a vacuum Bianchi II orbit in A 
which intersects the boundary of Ns in U. But then the limit sets of the orbit must also be in A. So 
for any neighborhood, Ne, of Ti, choose Ns with the following properties: 

i. Ng is a neighborhood of Ti in U. 

ii. No solutions in Ii are entirely contained in Ng and the unstable manifold of Ti is not entirely 
contained in Ng. 

iii. Ng C Ni, where iVi is a neighborhood of Ti which is topologically equivalent to Oi , as given at 
the end of section 4. 

iv. The limit points of all vacuum Bianchi II orbits which intersect Ng are contained in N^. 

None of these limit points is Ti. All of these limit points are non-flat Kasner points. At least one of 
these limit points is in A. □ 

6 Strong cosmic censorship 

While the results of this paper do not give a complete characterization of the a-limit set of each 
magnetic Bianchi VIq solution, they are sufficient for showing that none of these solutions has an 
extension past the initial singularity. That this is the case follows from the unboundedness of the 
Kretschmann scalar, Ra/s-ygR"^'*^ , in a neighborhood of the singularity. That the Kretschmann scalar 
is unbounded in a neighborhood of the singularity follows from the unboundedness of the expansion 
scalar in a neighborhood of the singularity and the presence of a non-flat Kasner point in the a-limit 
set (so there is a sequence of points on the orbit of any magnetic Bianchi VIq which both "goes to the 
singularity" and converges to a non-flat Kasner point). The argument is the same as in the proof of 
Theorem 5.1 in See also equations (3.5) and (3.6)]. The ratio between the Kretschmann scalar 
and the fourth power of the expansion scalar can be expressed as a polynomial in S+, A'^-i-, 
and H. This ratio is nonvanishing at a non-flat Kasner point, and the result follows. 

The same argument shows that the Kretschmann scalar is unbounded in a neighborhood of the 
singularity in each of the solutions considered in this paper which has a non-flat Kasner point in its 
a-limit set, and so such a solution cannot be extended past the singularity. This is the case for each 
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vacuum Bianchi VIq solution and each of the magnetic Bianchi 11 solutions considered in this paper. 
It is also the case for "most" of the vacuum Bianchi II solutions, the magnetic Bianchi 1 solutions and 
the Kasners themselves. 

The exceptions are the flat Kasner spacetimes, the rotationally symmetric vacuum Bianchi II 
spacetimes and the rotationally symmetric magnetic Bianchi I spacetimes. The a-limit set of each of 
these consists of a single flat Kasner point and the argument just given is inconclusive in these cases 
(the ratio between the Kretschmann scalar and the fourth power of the expansion scalar vanishes). 
It turns out that each of these spacetimes can, in fact, be extended past the initial singularity. It 
is already known that the flat Kasner spacetimes and the rotationally symmetric vacuum Bianchi II 
spacetimes can be extended. An extension for the rotationally symmetric magnetic Bianchi 1 
spacetimes is given in the appendix. 

7 Conclusion 

It has been shown in this paper that any spatially homogeneous solution to Einstein's equation of 
Bianchi type VIq with source consisting of a magnetic field orthogonal to the two commuting Killing 
vector fields has an oscillatory approach to the initial singularity and has no extension past the 
singularity. In the formulation which has been used, other spatially homogeneous solutions lie on 
the boundary of magnetic Bianchi VIq. All of the boundary solutions are convergent rather than 
oscillatory, and most are inextendible past the initial singularity. The only ones which are extendible 
are the rotationally symmetric vacuum Bianchi II solutions, the rotationally symmetric magnetic 
Bianchi 1 solutions and the flat Kasners. Since these are non-generic in the class, strong cosmic 
censorship is supported. 

The methods used here to characterize the approach to the singularity are those of dynamical 
systems analysis. One would like to obtain a rigorous proof that generic non-polarized magnetic 
Bianchi VIo solutions arc mixmaster. This result could be obtained by finding two functions which 
are strictly monotonic on all or at least generic magnetic Bianchi VIq solutions as formulated here 
and also, in turn, on the two classes of magnetic Bianchi type 11 solutions which lie on the boundary. 
The existence of such functions would exclude these solutions of the boundary from the a-limit sets of 
magnetic Bianchi VIq solutions. Alternatively, one could try to show that this class is mixmaster by 
other methods. It appears that such a result has been recently obtained by Ringstrom for solutions 
of vacuum Bianchi types Vlll and IX ||l^ . 

It is hoped that this work will contribute to the understanding of the nature of cosmological 
singularities in classical General Relativity. In particular, the validity in the generic case of the 
description of the approach to the cosmological singularity given by Belinskii, Khalatnikov and Lifshitz 
is still an open question. Whether or not the methods of dynamical systems analysis turn out to be 
directly useful in studying this question in spatially inhomogcneous classes of cosmological solutions to 
Einstein's equation, they have added to our understanding of the spatially homogeneous case, which 
gives a surer foundation for further study. 
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Appendix 



The metric of a rotationally symmetric Bianchi I spacetime with a magnetic field orthogonal to the 
plane of rotational symmetry can be written as follows Appendix C, Case 6]. 

g = -A^ df + t^A-^ dx^ + A^ dy^ + A^ dz^ (8) 

with A = 1 + jB^t^ . The Maxwell tensor is F = B {dy ® dz — dz ® dy), with B a constant. In the 
direction of the singularity (t 0) the generalized Kasner exponents converge to (1, 0, 0). This is the 
point Ti in the formulation used in the present paper. In the opposite time direction (t —^ oo) the 
generalized Kasner exponents converge to (— ^, |, |). This is the point Qi in the formulation used 
in the present paper. These solutions were given by Rosen |20) and studied by Jacobs but the 
author is not aware of a published extension of these spacetimes. 

Here is a coordinate transformation that allows each of these spacetimes (inequivalent spacetimes 
result from different values of the constant \B\) to be extended. 

i=t cosh X i = tsinha; (9) 

In the original spacetime, = — . But if a function D is defined hy D = — SP' , and if the 
function A is now defined A = 1 + jB^D, then the transformed metric 

g^A^ {-dP + d5p) + (-B^ + B^D^^—^^){S:di-idS:f + A^ dy^ + A^ dz^ (10) 

is a solution to Einstein's equation in the region A > 0. This is an extension of the original solution. 
The region defined by > £> > — -g^- is not in the original spacetime. The singularity at A = is a 
curvature singularity. Both the Kretschmann scalar and RapR"^ are unbounded as A — > 0. 

RcfS^fsR''^^' = (20 - 6B^D + \b^D^)^ R^pR'^P = ^ (11) 

di is everywhere timelike, but dx is not everywhere spacelike. However, g^g^x — 9^^ = ^1- Thus the 
two-dimensional surface defined by constant y and z is everywhere timelike. 
The singularity is timelike. 

Along the x-axis the spacetime becomes singular a finite distance from the origin in either direction. 
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Figure 1: The Kasner circle. The arc between T2 and T3 containing Qi is /Ci. On /Ci, < 0, while 
A+ > and A_ > 0. The arc between T3 and Ti containing Q2 is IC2- On JC2, A+ < 0, while A_ > 
and Xh > 0. The arc between Ti and T2 containing Qs is /C3. On /C3, A_ < 0, while A^ > and 
A+ > 0. 
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